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ABSTRACT 

Rosenau [Phys. Rev. A, 40 (1989), pp. 7193-6] has recently proposed a regularized 
version of the Chapman-Enskog expansion of hydrodynamics. This regularized expansion 
resembles the usual Navier-Stokes viscosity terms at law wave-numbers, but unlike the latter, 
it has the advantage of being a bounded macroscopic approximation to the linearized collision 
operator. . - , , , 

In this paper we study the behavior of Rosenau regularization of the Chapman-Enskog 
expansion (R-C-E) in the context of scalar conservation laws. We show that this R-C-E 
model retains the essential properties of the usual viscosity approximation, e.g., existence 
of travelling waves, monotonicity, upper-Lipschitz continuity etc., and at the same time, it 
sharpens the standard viscous shock layers. We prove that the regularized R-C-E approx- 
imation converges to the underlying inviscid entropy solution as its mean-free-path e J, 0, 
and we estimate the convergence rate_ r 
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1. Introduction 


Rosenau [R] has recently proposed the scalar equation 

u t + /(«) . = y: ~ ek 

(i.i) 


1 + m 2 e 2 k 2 
1 


“(*) 


= £ 


r*(*) 

J xx 


.1 + m 2 e 2 k 2 

as a model for his regularized version of the Chapman-Enskog expansion for hydrodynamics. 
The operator on the right side looks like the usual viscosity term eu xx at low wave-numbers 
k, while for higher wave numbers it is intended to model a bounded approximation of a 
linearized collision operator, thereby avoiding the artificial instabilities that occur when the 
Chapman-Enskog expansion for such an operator is truncated after a finite number of terms 
[R]. 

In this article we will compare the behavior of solutions of Rosenau regularization for the 
Chapman-Enskog expansion, abbreviated hereafter as the R-C-E equation, (1.1), with those 
of the viscous conservation law 


( L2 ) u t + f(u) x = eu xx , 

towards which (1.1) tends as m — ► 0, and with the conservation law with absorption 

C 1 - 3 ) u t + f(u) x = — 

era 2 

Since the right side of (1.1) tends to that of (1.3) in the limit of large k } it is not surprising 
that the smoothness properties of solutions of the former resemble those of the latter. In 
particular, the R-C-E equation does not smooth out initial discontinuities, but as shown 
in §2, it does preserve the smoothness of initially smooth small initial data. On the other 
hand, the right side of (1.1) also resembles that of (1.2) in that both are second derivatives. 
Consequently, it is plausible that the regularized R-C-E equation (1.1), like the ordinary 
viscosity equation (1.2), should have travelling wave solutions connecting shock states of the 
underlying conservation law 

C 1 - 4 ) u t + f(u) x = 0 . 

In §3 we show that when f n > 0 such solutions exist iff m is sufficiently small. 

At the same time, solutions of the R-C-E equation (1.1) also resemble those of the 
inviscid conservation law (1.4) in that both admit unique entropy solutions which share 
similar properties. In §4 we show that the R-C-E solution operator associated with (1.1), 
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like the entropy solution operator of (1.4), is Z^-contractive, monotone, and BV-bounded. 
Furthermore, the R-C-E solution of (1.1) tends to the inviscid entropy solution of (1.4) as 
the ’mean-free-path’ e j. 0. Finally, if f" > 0, the R-C-E entropy solution of (1.1) is also 
upper-Lipschitz continuous, in agreement with Oleinik’s E-condition which characterizes the 
entropy solution of (1.4), and in §5 we estimate the convergence rate of the former to the 
latter as e j 0. 


2. Smoothness 

It is well known that solutions of (1.2) are smooth for t > Oj i.e. initial discontinuities are 
smoothed out at positive times. In contrast, by looking at piece-wise constant initial data 
or at the linear case f(u ) = u, one sees that initial discontinuities of solutions of (1.3) are 
merely attenuated, not smoothed out, at positive times. Since the damping of (1.1) is less 
than that of (1.3), it is clear that (1.1) also does not smooth out initial discontinuities. On 
the other hand, if the (e-independent) initial data for (1.3) is smooth then it will remain so 
provided that m is sufficiently small (see below). The next Theorem tells us that the same 
holds for the R-C-E equation (1.1). 

Theorem 2.1. The solution of the R-C-E equation (1.1) remains as smooth as its initial 
data, 

(2.1) u(x,0) = «o(*), 
provided the initial data u 0 are sufficiently small so that 

(2.2) 2{m||uo||i,oo||/ ,, ('Uo)||L“} 1//2 + m2e ll/ W ( u °)IU“ll ti olk“ < 

Remark. Since we can ensure that (2.2) will be satisfied for any fixed initial data by 
making m sufficiently small, Theorem 2.1 can also be viewed as showing how the smoothness 
properties of the R-C-E equation (1.1) approach those of the viscosity equation (1.2) as 

m — ► 0. 

Proof. We will show formally that (2.2) implies a bound on the L°° norms of u and u*. 
Estimates for higher derivatives then follow in standard fashion; see [M]. Furthermore, this 
fact ensures that the formal estimates can be justified either by smoothing the initial data 
or by applying a further vanishing viscosity regularization. 
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The first step towards obtaining the desired bounds is to note that the right side of (1.1) 
can be written as 


(2.3) 


-ek 2 1 v 

-u(k) 


1 + m 2 e 2 k 2 




-1 f 

_ m 2 c \ U ~Q mf * U J » 
me 


where * denotes convolution, and 

(2.4) 
satisfies 

(2.5) 


<M*) = Y e e ~ w ' 


WQAW = 1 • 


To obtain a uniform bound on u, multiply (1.1) by | u p 2 |u and integrate over s; since 
\u p ~ 2 \uf(u) x is an exact derivative its integral vanishes, while the contribution of the right 
hand-side (2.3) is nonpositive, for by (2.5), 

' ,p- 2\ _i 


( 2 . 6 ) 


/ 7^r u ( u ~ Qme * u ^ dx ~ IHIL> 1 ll<3m e *'u|| L p} 


< 


-l 


m 2 e 


— IIQme L 1 } — 0. 


Dividing the remaining inequality by (p — 1 )||u|| L p and integrating over t from 0 to T, we 
obtain 


(2.7) 


Wr)IU> < ll“olltp, 


and the boundedness of ll«(r)|| Lf°° follows by letting p | oo. 

In order to estimate in similar fashion the L°° norm of u x , we differentiate (1.1), obtaining 


( 2 . 8 ) 


-1 


u xt + f(u) xx = — 5 -{u a - Q' m€ * u} , 


m £ e 


and as before, we multiply (2.8) by |u x | p 2 u x and integrate over x. Integrating by parts 
where necessary in the term containing /, and noting that by (2.4) 


(2.9) 


IIQJJU* = — , 


me 


we obtain after factoring out (p — l)||u x ||^ P , 


( 2 . 10 ) 


d.. .. 

j t \WAW + 


m 2 e 


m 3 e 2 


u & < 


{1 - m 2 £||/"(u)||L“(t,x)ll U x||L<»(t 1 x)}||'Ux||L y 
1 


m 3 e 2 


Ml 


Lr 
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Next we denote 


(2.11) Y(T) = m 2 e||/"('u 0 )||z,“’||«x(r)||i,«>. 

Applying Gronwall’s lemma to (2.10) and letting p | oo, we obtain that Y = sup t Y(T) 
does not exceed 

Y{r ) < Y < e [-T ^P 1 F(0) + yiy{l - )] }||«o||l« ||/ W ( u o)||l«> 

(2.12) 

< F(0) + Y^r||‘Uo|U~||/ ,, («o)IU“i 

as long as Y < 1. Estimate (2.12) is a quadratic inequality for Y for which the roots of the 
corresponding equation are 

(2.13) Y = -{1 + y(0) ± y/(l - ^(0)) 2 - 4m||txo||x,oo ||/ ,, (' u o) || l°° } • 

z 

Our assumption (2.2) tells us that the expression under the square root on the right is 
positive. Since y(0) is bounded by the smaller root in (2.13), it follows from (2.12) plus 
the continuity of Y(T) that Y(T) remains bounded by this root. This in turn confirms that 
indeed Y < 1, and the uniform bound of ||u*(T)||l« follows. □ 

Remark. Arguing along the above lines for the conservation law with absorption (1.3), 
one arrives at the inequality, analogous to (2.12), Y(T ) < e -T(1-r)/m e r(0). This shows 
that if y(0) < 1 then Y(T) remains < 1. Consequently, if Y^O) < 1, then Y(T) and hence 
iKcniii,-, satisfy a maximum principle in this case. 

3. Shock Profiles 

Lax’s generalized entropy conditions [L] for ’’legitimate” shock- wave solutions of the 
conservation law (1.4) can be interpreted as the requirement that these shocks can be realized 
as the limit of travelling wave solutions of the viscosity equation (1-2). If the flux function 
f is convex, these conditions reduce to the shock inequalities [L] 

(3.1) /'(«-) > 3 >/'(«+) > 

where s is the speed of the shock joining u_ on the left to u + on the right. In this section 
we show the analogous result for the (convex) R-C-E equation (1.1): It admits travelling 
wave solutions whose limit as e j 0 are shock wave solutions of (1.4), iff (3.1) holds and m 
is sufficiently small. 

Theorem 3.1. Assume f" > 0. Then (3.1) and the Rankine-Hugoniot shock condition 

(3.2) H(u + ) = 0, H(u) = -s{ u - u_} + {/(«) - 
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are necessary conditions for the existence of a travelling wave solution 

u(z = — - — ), lim,^ ±00 u(z) = u±, 

for (1.1). 

Conversely, if (3.1) ,(3. 2) hold, then a sufficient condition on m for the existence of such a 
travelling wave is 

(3.3) 4m 2 sup {— f"(u)H(u)} < 1 , 
and a necessary condition is 

(3.4) 4 m 2 {-f"(u.)H(u.)} < 1 . 

Here u, is defined by 


(3.5) 


f\uf) = s . 


Proof. Define z = and let ' denote Using (2.4) we find that a solution of (1.1) of 
the form u = u{z) satisfies 

(3.6) -,v! + f(u)' = {Q„ » u}" , 

where the convolution is now taken w.r.t. the variable z. The condition lim^iooU — » u± 
implies that also Q m * u — > u± as z — * ±oo, so there exist a sequence of values zf tending 
to ±oo on which ( Q m * u)' tends to zero. Hence, integrating (3.6) from zf to z and letting 
j —* oo we obtain 

(3.7) H{u) = -s{u-u_} + {f(u) - f(u_)} = {Q m *u}' . 

Now letting z tend to +oo cilong the sequence we find from (3.7) that H(u + ) = 0, i.e. 
(3.2) holds. 

Noting that H" = f" > 0, we see that H(u_) = 0 = H(u + ) implies H' = f — s < 0 at 
the lessor of u±, and H' > 0 at the greater of the two. Hence, if u + < u_ then (3.1) holds, 
while if this inequality is reversed then so are those of (3.1), i.e., we can replace (3.1) by the 
condition 

(3.8) u_ > u + . 
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Next, we apply to (3.7) the operator 1 — 771,2 (the inverse of the operator of convolution 
with Qm), to obtain 

(3.9) u' = {1 - m 2 -^}H(u) = H(u ) - m 2 {H'(u)u" + H"(u)(u') 2 } . 

We note that since all nonzero solutions g of {1 — m 2 gp-}</ = 0 are unbounded on i2, the 
solution of (3.9) with bounded u and v! which we construct below, also satisfies (3.7). To 
construct such a solution we introduce the auxiliary variable 

(3.10) v = u' , 
which enables us to rewrite (3.9) as the 2x2 system 

(3.11) u' = v 

(3.12) m 2 H'(u)v' = H(u) — v — m 2 H"(u)v 2 . 

The convexity of H{u) together with the Rankine-Hugoniot condition (3.2) imply that the 
only critical points of system (3.11-12) are (u_,0) and (u+,0). . 

We remark that the linearization of (3.11-12) near the critical points (u_,0) and (u+,0) 
shows that they are both saddles, so that topological methods (see e.g. [S] ) cannot be 
applied; one might even tempted to conclude that the existence of a trajectory joining these 
saddle points is unlikely. What saves the day, however, is the fact that the system is singular 
on the line u = u ,, i.e., that the coefficient H'{u) on the left of (3.12) vanishes at u*, which 
by (3.1) lies between u_ and u+. 

The key to finding a trajectory joining the two critical points is to note that solutions of 
(3.11-12) can cross the line u = u, only at points (u,,u*) which make the right side of 
(3.12) vanish: Equation (3.11) implies that H'(u(z )) is 0(z - z«) near the value z, for which 
u (z,) = u,, and hence (3.12) shows that |vj — » oo as z — ► z,, unless the right side of (3.12) 
tends to zero. Also, since the right side of (3.12) is quadratic in v, a comparison of (3.12) 
with the equations zv' = ±v 2 shows that in fact jv| reaches infinity before u reaches u*. 

In order to obtain a trajectory joining u_ at z = — oo to u + at z — +oo it is therefore 
necessary and sufficient to find trajectories joining (it_,0) (respectively, {u+, 0)) at z = — oo 
(respectively, z = +oo) to (u t ,v,) (with the same value of u, for both cases) at some finite 
values of z; we can always arrange for the two values of z to coincide because the system is 
autonomous. Since trajectories through (u*,u*) are not unique, the existence of our desired 
trajectories, which when put together join u_ to u + , no longer seems so unlikely. 
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Now the right side of (3.12) is a quadratic expression in v, whose roots are 


(3.13) 


v±(u) 


-l±y/l + 4m 2 ff"(u)i?(u) 

2m 2 H"(u) 


If the argument under the square root is negative at u = u, then clearly no such v, exists; 
this gives the necessity of (3.4) for the existence of travelling wave solutions. We now turn 
to discuss the sufficiency of condition (3.3): it says that v± exist for all u between u_ and 
u+, and we want to show that when this happens, then the trajectories mentioned above 
exist iff (3.8) holds. Namely, that these trajectories exist if we replace u_ and u + by 


(3.14) 


u =max{u_,u+} and u + = min{u_,u + } , 


but not with these latter two interchanged. 

The linearized system around the two critical points has the form 


(3.15) 



Since the determinant of the matrix on the right of (3.15) is negative, both critical points 
are saddles, as claimed. Now, it is not hard to calculate directly the asymptotic directions 
of the solutions that approach each critical point as z tends to ±oo, as these are simply 
the eigenvectors of the matrix in (3.15), but in any case we will have to determine from 
(2.11-12) the signs of v! and v' in various regions, and this information suffices to determine 
which region the various asymptotic directions lie in. In this way we obtain Diagram 1, the 
phase-plane diagram for the case when m satisfies (3.3). 

Based on Diagram 1, the existence of a travelling wave solution is argued as follows. There 
is a trajectory that leaves the critical point (u - , 0) and enters region I. If this trajectory 
remains in region I until u reaches the value u t) then by the above analysis it reaches the 
point (u t , v,); in this case, u' = v as well as u are monotonic on this semi-trajectory. The only 
way that the trajectory can leave region I before reaching the line u = u*, is by entering 
region II ; but v' > 0 in this region, so “clearly” the trajectory still reaches (u,,v*). A 
similar analysis backwards in the “time” z shows that there is a semi-trajectory from (u t ,v t ) 
to (tt + , 0). By checking the other trajectories leaving and entering each critical point we see 
that no trajectory joins u + to u~ or either point to itself. 

Although the above argument is sound provided that Diagram 1 is accurate, we have yet 
to verify one crucial feature of that diagram. Specifically, the argument assumed that if we 
enter region II at the point P and travel within this region keeping v' > 0 and u' < 0, then 
we cannot reach region III. (Clearly, no problem arises from the possibility of re-entering 
region I.). Thus we have to show that the situation shown in Diagram 2 is impossible. 
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Analytically, we must show that 

v-(ui) 5; u+(u 2 ) for u. < u\ < u 2 .legno(3.16) 

Defining Hi - - H{u { ) and A* = -4 m 2 H"(ui)H(ui), then (3.16) boils down to 

(3.17) —2Hi{l + s/l^A 1 }/A 1 < -2H 2 (1 - yfl ^A 2 }/A 2 . 

Now, the convexity of / (and hence of H ) together with the fact that H vanishes at u± imply 
that Hi > H 2 > 0; these facts together with our assumption (3.3) imply that 

(3.18) 1 > Ai > 0 . 

Therefore, a sufficient condition for (3.17) to hold is that 

(3.19) {1 + \/l — ~A\)IAx > {1 — ^/l — A 2 }/A 2 

for all Ai satisfying (3.18). A little algebraic manipulation shows that this is indeed the case. 
Consequently, (3.16) holds, i.e., no point such as the point Q in Diagram 2 can exist, and so 
the argument based on Diagram 1 is valid. 

As to increases past the value that makes equality hold in (3.3), we obtain the situation 
shown in Diagram 3. 

Namely, a gap appears in region II, through which our trajectory might possibly plunge 
into the abyss of region III. Hence we cannot say whether a trajectory joining u~ to u + 
exists or not. Finally, when m increases past the value that makes equality hold in (3.4), 
then the phase-plane looks like Diagram 4, and the descent of our trajectory to — oo becomes 
a certainty. □ 

We close this section by quantifying Rosenau’s statement [R] that the travelling-wave 
solutions of the R-C-E equation (1.1) give narrower shock layers than those of the viscosity 
equation (1.4). We will adopt as our measure of shock width, w = (u_ - u + )j\u ' t | with u' 0 
evaluated at the point u, at which H'(u) = 0. (It should be noted, however, that although 
this value of u[ is always maximum for (1.4), it will be maximum for (1.1) only if the 
trajectory does not enter region II of Diagram 1, which in turn is guaranteed by our analysis 
only when the curve v — u + (u) has its unique local minimum at tt*.) Since the relevant value 
of u ^ for (1-4) is given by v! n — — H{^ u,), while the value of u ^ for (1.1) is ' U +('U , *)> the estimate 

/g go) I < ^Chapman-En.kog < j 

2 U> v ; 5C0ua 

follows from the simple lemma; If a quadratic equation has real roots, then the root closer 
to zero lies between r and 2r, where r is the root of its linear part. 
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4. Entropy Solutions and the Zero Mean-Free-Path Limit 

The parameter m does not play a role in our analysis in this section, and so will be set 
equal to 1 for convenience. 

Since solutions of the R-C-E equation (1.1) may contain singularities weak solutions must 
be admitted. Since the latter need not be unique, we single out an “entropy” solution of the 
R-C-E equation (1.1) as the one satisfying the Krushkov-like [K] inequality 

d t \u e -c\ + d x [sgn[u e - c){f(u e ) - /(c)}] 

(4.1) 

< ~j{\u e - c\ - sgn(u e - c)Q e * (u e - c)} , 

for all real c’s. In particular, by choosing c = +sup|u,| (respectively, c = — sup |u*|), we 
obtain from (4.1) that u e is a supersolution (respectively, a subsolution) of (1.1), and hence 

(1.1) is satisfied in the sense of distributions. We turn to show that (4.1) admits a unique 
solution u t , and that this solution converges to the unique entropy solution of (1.4) as e goes 
to zero. 

Theorem 4.1. For any uo in BV there exists a unique solution u e of the R-C-E equa- 
tion (4-1) i and as £ j 0, u e converges in L 1 to the unique entropy solution of the inviscid 
conservation law (1-4)- 

Proof. 

Add the artificial viscosity term 6u xx to the right side of (1.1); the resulting equation has 
a unique smooth solution u s e . By a straightforward adaptation of Krushkov’s [K, section 4] 
proof for the artificial viscosity method for (1.4), we obtain that the set {u s e }s> o is bounded 
in BV (uniformly in £ and £) and precompact in L 1 , and hence that a subsequence converges 
as 8 — > 0 to a solution u e of (4.1). Similarly, by the argument on pages 224-5 of [K] we 
obtain from (4.1) the consequence 

/ / {|u e - v*|$, + sgn(u e - v e )[f(u e ) - f(v t )\$ x }dxdt 

JO J — oo 

( 4 - 2 ) 

pT POO 1 

> / -{|u e - v e \ - sgn(u e - v e )Q t * (u e - v e )}$dxdt 

Jo J— oo € 

where $ is an arbitrary nonnegative test function. 

Next, we remark that the expression inside the curly brackets on the right of (4.2) need 
not be positive, but in view of (2.5), its spatial integral is. Therefore, by choosing $(<, x) = 
$i(i)$ 2 (z), anc ^ letting $2 tend to the function that is identically one, we obtain 

(4.3) [ f \u € - v e \$ u > 0 . 

Jo J — oo 
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Continuing as in [K], we let approach the indicator function of the interval [0, t] to 
conclude 

(4.4) ||ti«(t) - v e (t ) || L i < ||u,(0) - Ve (0)|Ux . 

In particular, this shows that the solution of (4.1) is unique. 

The solutions {it,} of (4.1) inherit the BV bound of the {u*}, and the argument of section 
4 of [K] shows that this bound implies precompactness in L 1 . Hence as e — ► 0 a subsequence 
converges to a weak solution u of (1.4). Because the right side of (4.2) is known to be positive 
only when $ has no dependence on x, we cannot use the entropies of (1.4) (as in [Ta]) to 
conclude that u is the entropy solution of (1.4). However, (4.4) implies the corresponding 
estimate for the weak solutions u and v obtained in the limit as e goes to zero, and by an 
argument of Lax [L] this suffices to show that we obtain the entropy solution: It is not hard 
to see that when (1-4) has a smooth solution then our scheme must converge to that solution. 
Hence by the corollary to theorem (3.5) of [L], any solution u of (1.4) obtained in the limit 
e — ► 0 from (4.1) has the property that all of its discontinuities satisfy the generalized Lax 
shock inequalities. By theorem (3.5) of [L], this implies that u is the unique entropy solution 
of (1.4). Finally, since any sequence of cs tending to zero has a convergent subsequence, the 
uniqueness of the limit shows that convergence holds without passing to a sequence. □ 


5. The Convergence Rate of the Zero Mean-Free-Path Limit 


Theorem 4.1 shows that the R-C-E equation (1.1) retains several properties of the viscous 
conservation law (1.2). In particular, (4.4) asserts that the solution operator is an L 1 - 
contraction, and hence by conservation plus translation invariance it is monotone [CM, 
Lemma 3.2], and by translation invariance it is BV-bounded: 


(5.1) 


K,(<)||bV < ||“e(0)||flV- 


Next we show that the nonlinear R-C-E equation (1.1) also satisfies Oleinik’s E-entropy 
condition, e.g., [Sm],[T]. 

THEOREM 5.1. Assume f" > a > 0. Then the following a priori estimate holds* 

1 


(5.2) 


u e (t)IUi P + < 


U e(°)lll„>+ + at 


t> 0. 


’We let H^llLtp, IMUi P + and ||^||x,i p ' denote respectively esssup,^ | j , esssup I9 £ V [^r^] + * 
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Remark. The inequality (5.2) implies that the positive-variation and hence the total- 
variation of u e (i) decays in time. Furthermore, this gives us another proof of the zero 
mean-free-path convergence to the entropy solution of (1.4) for any Lj^-initial data uo (cf. 
Corollary 5.2). 

Proof. We add the artificial viscosity term Su xx to regularize (1.1), obtaining 

(5.3) d t u i + d x f(u s e ) = -Ar e {< - Q m . * <} + Sdlu 6 ' . 

Differentiation of (5.3) yields for w = d x u 6 e , 

(5.4) d t w -1- f\u s )d x w + f"(u s )w 2 --= — ^-{to - Q me * + Sdlw. 

me 

Hence, since f > a > 0, it follows that W(t) = max x w{i) is governed by the differential 
inequality 

(5.5) W{t) + aW 2 (t) < -Qme*W}< 0 

me 

and (5.2) follows by letting 5 [ ().□ 

Theorem 5.1 shows that solutions of the R-C-E equation (1.1) are Lip + - stable. Moreover, 
(5.1) implies that the Lip 1 - size of their truncation if of order 0(e), for 

(5.6) \\d t u e + d x f{u e )\\ Up ' = e\\Qmc * d x u e \\ L i < e||<? me || L i|Mt)||flv < £||u.(0)||bv. 

Using the result of [T] we conclude that the Lip'- convergence rate of the R-C-E solutions 
to the corresponding entropy solution is also of order 0(e). 

Corollary 5.2. Assume that f" > a > 0, and let u e be the unique R-C-E solution of (4-1) 
subject to C 1 initial conditions u e ( 0) = / . then u c converges to the unique entropy solution 

of (1.4) and the following error estimates hold 

||u e (t) — ■u(t)||w'-'.J> < Const ■ £ , 1 < p < OO, .8 = 0,1. 

Remark. The choice (s, p) = (1, 1) corresponds to Lip 1 -convergence rate of order 0(e). The 
choice (s,p) = (0, 1) corresponds to the usual viscous i 1 -convergence rate of order 0(e*). 
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